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QUASIPARTICLE  DAMPING  OF  A POSITRON  IN  AN  ELECTRON  GAS 

J . Oliva 

Department  of  Physics,  University  of  California  San  Diego 
La  Jolla,  California  92093 


ABSTRACT 

We  present  a calculation  of  the  imaginary  part  of  the  self  energy  £j 
of  a positron  in  an  electron  gas.  £j  is  inversely  proportional  to  the 
inelastic  collision  rate  between  the  positron  and  the  electron  gas. 

The  calculation  was  done  using  standard  many  body  techniques  and  within 
the  RPA  approximation.  There  are  two  contributions  to  £ : one  arising 
from  particle-hole  pair  excitation,  the  other  from  plasmon  excitation. 
The  former  dominates  at  low  energy  where  - E2  while  the  latter  dooi- 


ABSTRACT 

We  present  a calculation  of  the  imaginary  part  of  the  self  energy  Ej 

of  a positron  in  an  electron  gas.  Ej  is  inversely  proportional  to  the 

inelastic  collision  rate  between  the  positron  and  the  electron  gas. 

The  calculation  was  done  using  standard  many  body  techniques  and  within 

the  RPA  approximation.  There  are  two  contributions  to  Ej : one  arising 

from  particle-hole  pair  excitation,  the  other  from  plasmon  excitation. 

The  former  dominates  at  low  energy  where  Ej  ~ E^while  the  latter  domi- 

1 E 

nates  at  high  energy  where  Ej ^ £n  — . 


INTRODUCTION 


The  collision  rate  of  a positron  moving  in  an  electron  gas  is  investigated 
using  a self-energy  approach.  The  quasi-particle  self-energy  is  a complex- 
valued  function  of  momentum.  The  imaginary  part  is  a manifestation  of  real 
collisions  between  the  positron  and  the  electron  gas.  In  Section  I,  a Green's 
function  for  a positron  in  an  electron  gas  is  introduced.  The  self  energy  £ 
is  identified  and  the  quasiparticle  collision  rate  is  then  given  by  an  approxi- 
mate expression  involving  its  imaginary  part  £j.  In  Section  II  an  examination 
of  £j  within  RPA  is  carried  out.  Two  contributions  to  £j  are  found  and  inves- 
tigated; one  arising  from  particle-hole  pair  excitation  (Section  III),  the  other 
arising  from  plasmon  excitation  (Section  IV). 


I.  POSITRON  GREEN'S  FUNCTION 


We  consider  a uniform  system  of  a large  number,  N,  of  electrons  and  a 
single  positron,  all  moving  in  a uniform  neutralizing  background  and  interacting 
via  the  Coulomb  force.  The  second  quantized  Hamiltonian  for  the  system  may  be 
written  in  the  Schrodinger  picture  as* 


H * H 
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* 1 throughout  the  discussion. 
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where  ^(x)  0J>a (x) ) are  Schrodinger  electron  field  operators  which  create 

(destroy)  an  electron  of  spin  a at  the  position  x,  and  where  <b+(x)  (6  (x))  are 

~ a - a ~ 

Schrodinger  positron  field  operators  which  create  (destroy)  a positron  of  spin 
a at  the  position  x.  m denotes  the  mass  of  the  electron  and  positron. 

We  seek  to  calculate  the  lifetime  of  the  positron  quasiparticle  states. 

To  this  end  we  introduce  a positron  Green's  function  GPg(xt,x't')  appropriate 
to  the  description  of  propogation  of  a single  positron  in  an  electron  gas  and 
defined  as  follows: 

-1  AiI<W;t)-4<i,t'))i',,o> 


GP  (xt.x’t*) 
ap  - - 


(2) 


P 


where  <j>  (xt)  and  <J>uD(x't')  are  the  Heisenberg  positron  field  operators 
and  where  j denotes  the  Heisenberg  ground  state  of  a uniform,  fully 

° M 

interacting,  N-electron  gas.  T denotes  the  Dyson  chronological  operator  which 
orders  operators  in  accordance  with  decreasing  time  arguments  toward  the  right. 
Since,  in  our  case,  H does  not  involve  the  spin  variables,  GPg  is  diagonal  in 
the  spin  indices;  we  only  speak  of  the  diagonal  part  GP(xt,x't)  from  here  on 
out  and  supress  spin  indices  where  appropriate. 

That  the  Green's  function  as  defined  in  Eq.  (2)  is  pertinent,  is  suggested 
in  the  corresponding  Lehmann  representation.  To  obtain  this  we  introduce  a 
completeness  relation  in  Eq.  (2) 

GP(xt,x'f)  - -iH  |9(t-t,)('yo|<t>H(xt)|'i'n)('Fnl<t>g(x,t,)|'l'o> 

- e (t'-t)(rj  **(;’  t ' ) | 'QCi' J *H(;t>  I 

(3) 
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where  the  intermediate  states  are  eigenstates  of  the  Hamiltonian  H | Eq  - (1) 

for  the  fully  interacting  hybrid  system  of  an  arbitrary  number  of  electrons  and 
either  zero  or  one  positron;  and  where,  as  such,  the  summation  in  Eq.  (3)  in- 
cludes states  of  all  possible  numbers  of  electrons  and  either  zero  or  one 
positron.  Note  we  have  taken 1* 

Since  contains  no  positrons,  the  second  term  of  the  summand  of  Eq.  (3) 

vanishes  for  all  n.  Moreover,  only  intermediate  states  with  N electrons  and  a 
single  positron  will  contribute  in  the  remaining  summand. 

Upon  invoking  the  spatial  and  temporal  translational  invariance  of  the 
system,  the  following  Lehmann  representation  for  the  Fourier  transform  of 
GP(xt,x't)  is  readily  obtained: 


Gp(k,aj)  = L'  \ p 
n ~ ’in 


w - (E  -E)  + in 
n 


(4) 


where  the  prime  on  the  summation  denotes  the  restriction  to  summation  over 

intermediate  states  with  N electrons  and  one  positron,  where  P and  E are  the 

r ~n  n 

momentum  and  energy  of  the  intermediate  state  | where  E is  the  energy  of 

the  fully  interacting  N-electron  ground  state  and  where  <J>(o)  and  <J>+(o)  are  now 

the  Schrodinger  positron  destruction  and  creation  field  operators,  evaluated  at 

the  origin  of  coordinate  space.  6.  p is  the  Kronecker  delta  function. 

- ’ ~n 

Note  importantly  that  the  poles  of  G(k,aj)  in  the  co-plane  occur  at  the 
exact  values  of  the  change  in  energy  associated  with  the  addition  to  an  electron 
gas  originally  in  the  ground  state  of  a single  positron  such  that  the  final 
system  has  a momentum  k.  Note,  by  way  of  example,  that  when  the  interaction 
between  the  positron  and  the  electron  gas  is  set  equal  to  zero,  the  intermediate 
states  | ¥ appearing  in  Eq.  (4)  separate  into  a product  of  a free  positron  state 
vector  times  an  N electron  (fully  interacting)  state  vector,  leading  to 
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(5) 


- 


G F(k,u)) 


k2 

2S  + irl 


k d 

being  the  energy  of  a free  positron.  This  form  for  G p is  obviously  consistent 

with  the  previous  statement  regarding  the  positions  of  the  poles. 

Now  we  would  like  to  describe  these  many  body  excitations  in  terms  of  a 
positron  quasi-particle  picture.  Quasi-particle  states  are  long  lived  quasi- 
eigenstates of  the  many  body  system  which  have  properties  approximating  those 
of  free  single  particle  states,  though  differing  from  these  in  being  damped 
with  time,  this  reflecting  the  presence  of  interaction  with  the  rest  of  the  many 
body  system. 

In  light  of  the  form  Eq.  (5),  we  are  led  in  the  standard  fashion  to  the 
following  quasi-particle  approximation  for  the  fully  interacting  GP: 


GP(k,w) 


Z(k) 


ZOO 


to  - E (k)  + in  to  - E(k)  - iT(k) 


(6) 


where  in  the  last  step  we  have  written  the  complex  valued  quasi-particle  energy 
E (k)  explicitly  in  terms  of  its  real  part,  E(k),  and  its  imaginary  part,  T(k). 
Z(k)  is  an  unimportant  function,  related  to  the  strength  of  the  quasi-particle 
state. 

The  quasi-eigenstates  behave  in  time  as  e Thus  the  probability 

density  of  the  quasi-particle  having  initial  excitation  energy  E(k)  will  decay 
in  time  exponentially  at  a rate  y given  by: 


y(k)  -a|roo| 


(7) 


This  can  be  interpreted  as  a total  collision  rate  of  the  positron  in  the  electron 
gas,  and  the  inverse  therefore  viewed  as  a lifetime.  We  point  out  that  for  the 
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•quasi— particle  concept  to  be  useful  the  damping  rate  must  be  small  compared 


to  the  excitation  energy: 


TCk)  « E (k) 


Having  established  the  relevance  of  GP,  we  may  now  bring  to  bear  on  the 
matter  of  determining  the  positron  quasiparticle  collision  rate  the  standard 
methods  of  modern  many  body  theory.  We  rewrite  Eq.  (2)  in  terms  of  the  inter- 
action picture  via  the  Gell-Mann-Low  theorem,  thereby  obtaining  the  following 
perturbation  series: 


iGp(xt,x't')  = ^o|T(4,R(xt)  <j)J(x't')|^o}  = <*0|S|*0>  1.  (*el£(-i) 


H I 


J dt /.  e ’ * 1%^ 


^f.|t(HINT^(ti) 


• • H!NTI(tV)  *i(?c) 


where  the  subscript  I on  state  vectors  and  operators  denotes  interaction  picture. 

is  the  K-electron  non-interacting  ground  state  and  HINT  (t^)  is  obtained 
from  Eq.  (lc) . S is  the  scattering  operator,  i.e. , S » U^O*,-00),  Uj(t.t')  being 
the  conventional  time  development  operator  (which  takes  the  form  of  the  numerator 
of  Eq.  (9)  with  the  two  positron  field  operators  absent). 

Using  the  explicit  form  of  in  terms  of  field  operators,  this  perturba- 

tion series  may  be  expressed  more  conveniently  by  applying  Wick's  theorem  to 
evaluate  the  expectation  values  of  the  time  ordered  products  of  field  operators. 


er.  1 


The  resulting  form  has  a diagrammatic  interpretation  completely  analogous  to 
that  of  the  more  familiar  case  of  the  single-particle  Green's  function  of  an 
electron  in  an  electron  gas,  and  with  completely  analogous  Feynman  rules.  As 
in  the  just  mentioned  case,  it  is  readily  shown  that  the  "disconnected"  diagrams 
cancel  the  denominator  in  Eq.  (9)  leaving  only  the  connected  diagrams. 

We  illustrate  a few  of  the  lower  order  diagrams  resulting  from  this  procedure: 


A double  (single)  simple  line  denotes  a free  positron  (electron)  Green's  function; 
a double  (single)  wavy  line  denotes  a bare  positron-electron  (electron-electron) 

4 

Coulomb  interaction. 

It  is  clear  that  we  may  arrive  at  a Dyson's  equation  for  GP(xt,x't)  which 
will  involve  a positron  proper  self  energy  EP(xt,x't'),  this  latter  being  the 
sum  of  all  irreducible  subdiagrams  arising  out  of  the  perturbation  expansion  for 
GP.  The  subdiagrams  referred  to  result  from  removing  the  "ends"  of  the  diagrams 
for  GP.  ^The  irreducible  subdiagrams  are  those  which  cannot  in  turn  be  separated 
into  two  subdiagrams  by  removal  of  a single  particle  line.j  Invoking  the  spatial 
and  temporal  homogeneity  of  the  system,  we  may  Fourier  transform  Dyson's  equation, 
arriving  at: 

o 0 * n 

Gp(k,u)  = G P(k,oi)  + G P(k,ai)  E (k,ui)  Gp(k,u>). 


(10a) 


GP(k,u))  = * 

u - E°(k)  - E P(k,w) 


(10b) 
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o k 

where  E (k)  is  the  free  particle  energy  -r—  . 

Zm 


Upon  comparing  Eq.  (6)  and  Eq.  (10b),  we  see  that  we  can  now  relate  the 
quasiparticle  energy  E (k)  to  the  proper  self  energy  E (k,u>)  |we  here  on  out 
delete  the  superscript  p (denoting  posi*--on)  from  G & 5lj  . in  the  case  of 
long-lived  excitations: 


E(k)  sr  E(k)  + Re  E (k,E(k)) 


(11a) 


r(k)  * (1  " Is  Re  E’'(k,w)  l(jJ=E(k)^~1  ^O^OO) 


(lib) 


Our  interest  here  is  in  the  damping  rate  | 2r | , We  will  make  the  approximation 
that  the  derivative  term  appearing  in  Eq.  (lib)  is  small  compared  to  unity, 
giving: 


T(k)  Im  E (k,E(k)) 


(12) 


This  is  certainly  valid  for  very  large  k,  since  E (k,E(k))  0 for  k °°, 

and  is  likely  to  be  valid  for  very  small  k,  where  Re  E (k,E(k))  should  be 
relatively  flat.  For  intermediate  k,  this  approximation  is  questionable  though; 
we  shall  return  to  this  point  later. 

II.  APPROXIMATE  EVALUATION  OF  IM  £* 

We  now  evaluate  the  positron  proper  self  energy  in  an  electron  gas  within 
the  random  phase  approximation  (RPA) . By  this  it  is  meant  that  the  proper 
self  energy  is  approximated  by  an  infinite  summation  of  bubble-type  diagrams 
of  the  form  shown  in  the  following  Eq. 


2T  - r 


>vw 


rv 

L 


* 
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(13) 


We  have  introduced  with  an  obvious  diagrammatic  identification  the  effective 

interaction,  of  a positron  in  an  electron  gas  . These  diagrams  for  Z (p,w) 

correspond  to  successive  virtual  excitations  of  electron  particle-hole  pairs 

by  the  positron.  Though  it  is  true,  as  with  the  case  of  electron  RPA,  that 

* 

this  approximation  for  Z becomes  very  good  (exact,  in  fact)  only  in  the  extreme 

high  density  limit  r « 1 Jr  being  the  Wigner-Seitz  radius  = mean  distance 

s i s 

between  electrons j,  we  probably  can  expect  qualitative,  even  semi-quantitative 
validity  for  rg  in  the  metallic  range  -2  <rg  <~5. 


We  note  that  the  exchange  diagram. 


, appearing  in  the  electron  Z 


has  no  analog  here  since  there  is  only  one  positron  in  the  system.  Furthermore, 
since  in  all  diagrams  for  there  are  exactly  two  electron-positron  bare 

Coulomb  interaction  lines,  each  term  is  equal  in  value  to  the  corresponding 
term  in  the  effective  potential,  V®^,  seen  by  an  electron  (e.g.,  3 = 3' 
etc).  We  thus  have 

vPff(k,u>)  - V*ff(k,u>)  - v (k) 


where  is  the  propagating  RPA  (Lindhard)  dielectric  function  of  an 

N-electron  gas  |defined  below,  Eq.  (28)j  and  where  v(k)  is  the  bare  Coulomb 


interaction: 
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v(k)  = 


4TTe 


(15) 


The  positron  proper  self  energy  is  given,  in  accordance  with  the  Feynman 
rulee,  by: 

Z*  (p,u>)  = — f d3k  du>'  G°(p-k,  uwd')  Vp  (k,u>’)  (16) 

RPA  (2it)4  J ~ ~ eff  ~ 

which,  using  Eq.  (5)  and  Eq.  (13),  becomes 


- — 


(2tt) 


J d3k  J du' 


v(k) 


[eKPA(k,a)'\  f w-(jJ,-E0(p-k)  + in 


(17) 

The  frequency  integral  appearing  in  Eq.  (16)  is  along  the  contour  C |shown  in 

Figure  (la)j,  so  chosen  so  as  to  properly  avoid  the  poles  of  £(k,o>'). 

RPA 

We  have  seen  that  in  our  approximation,  the  imaginary  part  of  the  self 
energy  is  directly  related  to  the  collision  rate  | Eq.  (12) j . A rather  convenient 

form  for  Im  now  obtained  we  from  now  on  omit  the  subscript  RPA  on 

* 1 

Z and  £RpAj-  Note  the  free-positron  G-F  in  the  integrand  of  Eq.  (17)  has  a 
RPA  J 

pole  at  to'  = oj  = a)  - E°(p-k)  + in.  Note  further  that  since 


(k,u>*) 


- 1 ~ 


0) 


1_ 

>2 


CO' 


(18) 


we  may,  by  Jordan's  lemma,  close  the  contour  C in  Eq.  (17)  with  the  standard 
infinite  semicircle  in,  e.g.,  the  upper  half  plane.  By  rotating  the  resulting 
closed  contour  90°  counterclockwise  |see  Fig.  (lc)J , we  have  by  Cauchy’s  theorem 
and  upon  noting  that  for  both  initial  and  final  contours,  the  integral  over  the 
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semicircular  part  vanishes. 


E (p  .a) ) 


d3k  du'  v(k) 


^ e (k,a/)  " 


(D-oi'-E  (p-k)  +in 


-?vi  / d3k  ( - 1)  v(k) ! 

0<E°(p-k)<w  \ (k,uj-E  (p-k)  +in)  / 

« w 'v  ** 


where  C'  is  the  +90  rotated  version  of  C (with  the  dogleg).  Realizing  that 
there  are  no  poles  of  the  first  integrand  in  Eq.  (18)  near  the  imaginary  axis, 
we  may  now  equivalently  replace  C'  by  the  imaginary  axis  and  in  addition,  set 
n*  0 in  the  integrand.  The  second  term  in  Eq.  (18)  is  a residue  contribution 
which  arises  since,  as  the  original  contour  is  rotated,  the  pole  u>^  of  the 
positron  propogator,  when  it  is  in  the  first  quadrant  |i.e.,  when  k^-u>  > E(p-k)  : 
note  limits  in  the  second  integral  of  Eq.  (1  )|  goes  from  being  within  to  being 
outside  the  closed  contour. 

It  can  be  shown  that  the  imaginary  part  of  the  line  integral  appearing  in 
Eq.  (1  ) vanishes  paving  replaced  C by  the  real  axis  and  setting  ttfO,  introduce 
the  change  in  variables  u*i(J.  Noting  then  that  since  e(k,<*)' ) is  the  propogating 
dielectric  function  which  has  the  symmetry  property  £(k,w')  = £(k,-w'),  the 
resulting  integral  is  manifestly  realj. 

Therefore  we  arrive  at  the  useful  expression 


Im  E* ( p,w)  = — Im  / d3k  

2 IT2  0<E(p-k)<u>  (k,u>-E0(p-k)  + in) 


which  will  serve  as  the  starting  point  for  the  subsequent  development. 


We  now  review  the  properties  of  the  Lindhard  dielectric  function  e(k,u>) 


Its  real  and  imaginary  parts  Er(k,oj)  and  E^k.ui)  are  given  by 


RPA 


where  pQ  is  the  Fermi  momentum  of  the  electron  gas. 
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n being  the  electron  density.  Note  we  have  introduced  a separate  notation 
for  the  two  distinct  analytic  forms  of  e^,  there  applying  in  different  regions 
of  the  u)  vrs.  k plane. 

We  illustrate  the  single-particle  excitation  spectrum  of  a free  electron 
gas  in  Figure  (2).  The  shaded  band  corresponds  to  the  region  where  excited 


particle-hole  pairs  can  exist.  Note  that  e^k.aj)  is  non  zero  only  within  this 


band.  The  two  bounding  curves 


,2  k p 

.w  - f + — 

max  2m  m 


(23a) 


, 2 k p 

w(k)  1° 

min  2m  m 


(23b) 


give  the  maximum  and  minimum  energies,  resp.,  of  a particle-hole  pair  of 
momentum  k. 

Also  shown  in  this  Figure  is  the  plasmon  excitation  branch  u)p(k)  which 
appears  when  interactions  are  allowed  for.  a)p(k)  solves  the  equation  of  the 
longitudinal  resonance  condition: 


e(k,  u)p(k) ) =*  0 


This  branch  starts  at  k=0  with  energy  uiq  given  by 


Air  n 

“o  -(-T-J  • 


rises  to  a maximum  value  (i)p(kc)  at  the  so-called  plasmon  cutoff  momentum,  kc> 

and  ceases  to  be  a well  defined  excitation  for  k > k . The  cutoff  is  determined 

c 

by  the  condition  that  it  becomes  possible  for  a plasmon  to  decay  into  a single 


pair  excitation,  i.e.,  when 


k2  p0 

u>  (k)  = -r—  + k — = (1)  (k) 

max  2m  m p 


T 


* 


Returning  now  to  the  calculation  for  Im  Z , it  will  be  seen  that  there  are 
in  general  two  types  of  contributions  to  Im  I*  jsee  Eq.  ( 20) | : one  from  the 
integration  over  the  region  of  k-space  where  e^(k,  oj  -E(p-k)  + in)  is  finite 
|this  corresponding  to  the  positron's  exciting  particle-hole  paris j and  the 
other  arising  from  any  poles  in  the  integrand  (i.e.,  zeros  of  falling  within 
the  domain  of  integration  Jthis  corresponding  to  the  positron's  exciting  plasmons 
Before  proceeding,  we  make  one  further  approximation:  we  take  for  the 
approximate  collision  rate  2r(see  Eq.  (13))  twice  the  value  of  the  self-energy 
E(k,E(k)),  evaluated  at  frequency  F(k)  • E°(k),  i.e.,  the  on  shell  value.  We 
will  thus  be  evaluating. 


Im 


Z*(p) 


2 TTZ 


Im  J 


d3k 


v(k) 


0<E°(p-k) <E°(p) 


£(k,  E°(p)-E°(p-k)  +in) 


(27) 


I 


III.  PARTICLE- HOLE  CONTRIBUTION 

As  indicated  in  Eq.  (21),  the  imaginary  part  of  the  dielectric  function 

1 2 I i 2| 

e(k,to)  is  non-zero  when  either  pQk  + k _>  mu  _>  I p^k  - ^ k I or 

1 2 

k <2pQ  & 0 < mu _<  kpQ  - ^ k . Using  the  fact  that  in  the  integral  of  Eq.  (27), 
e(k,ii>)  enters  with  frequency  argument  <i>=  (2pk  cos  9-k1 2)  jwhere  9 is  the 
angle  between  j>  and  k|,  it  is  straightforward  to  show  that  this  set  of  inequalities 
reduces  to  the  one  simple  inequality: 

|k  - p cos  0 1 _<  pQ  (28) 


The  region  of  k-space  specified  by  this  inequality  is  bounded  by  a surface 
known  as  a Pascal  limacon  of  revolution.  For  p > pQ,  there  are  two  branches  of 
the  surface  (one  contained  in  the  other)  whereas  for  p < pQ,  there  is  a single 
brar.ch  fsee  Fig.  (3)j  . The  region  of  k-space  in  which  Im  £(k  • ^ (2pk  cos9-k2)) 
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is  nonzero  and  which  contributes  to  the  integral  Eq.  (27)  is  then  the  region 
which  is  both  within  the  Pascal  limacon  (when  p > pQ  this  means  between  the  two 
branches)  and  within  the  sphere  of  radius  p centered  at  k * p.  |see  Fig.  (3)|  . 

•k  it 

We  find,  defining  E , (p)  = particle-hole  contribution  to  Im  E (p) , 

p-h 


Z.e 

2 


..  r r2px 

— Im  I dx  1 dk  — 

IT2  J J e(k,-r- 

o*o  r 2m 


(2pkx-k  ))+ei  (k.-^jjUpkx-k2)) 

(P<  PQ) 


(29a) 


ZT  n h(P> 
I,  p-h 


2 TT 


(2pkx-k2))+e1  (k,^(2Pkx-k2)) 


r 1 rPX+P 
dx 


_o 

P 


*o 

I dk  

I e (k, 

Jpx-p„  rv 


-^(2pkx-k2)+e  (k,^(2pkx-kzj) 


(P  >P0) 


(29b) 


We  first  consider  the  case  p « pQ.  Entering  into  Eq.  (29a)  is  then  the 
dielectric  function  evaluated  for  momentum  and  energy  very  small  compared  to 
PQ  and  Eq,  respectively: 


e(k,  "2^(2pk  cos9-k2))  ~ <1  + 


2mpc 

TT  k2 


+ i 


[2m2Po  /l 

I \2m 


— (2pk  cos0-k  ) 


(k  « pQ) 


(30) 


15 


see  Eq.  (21).  We  readily  arrive  at: 


Cr  * (P) 
p-h 


-ifer 


p <<  p. 


We  note  that  in  this  small  momentum  regime,  the  p-h  contribution  to  the 


collision  rate  decreases  with  increasing  electron  density.  Contributing  to 
this  is  the  fact  that  at  higher  densities,  the  characteristic  size  of  the 


region  of  enhanced  electron  density  screening  the  positron  decreases,  rendering 
the  positron  quasiparticle  less  effectively  scattered,  and  the  fact  that  at 
higher  densities,  where  the  kinetic  energy  dominates  the  dynamics,  the 


electrons  are  less  able  to  respond  to  the  presence  of  an  external  charge. 

In  passing,  we  note  the  obvious  physical  correspondence  between  an  electron 

of  momentum  kg  just  outside  the  Fermi  sphere  and  a positron  of  very  small 

momentum  k^  compared  to  pQ,  in  that  the  respective  quasiparticle  states  are 

both  very  long  lived  (t  - (k  -p  )2  for  an  electron.) 

e e o 

We  now  consider  the  case  of  high  momentum  p » pQ.  Starting  with  Eq.  (29b) 

* 

it  can  be  shown  that  for  this  range  of  p,  E (p)  varies  inversely  with  p: 

P-h 


C-h(F)  * f C(Po> 


(32a) 


C(pQ)  = Im 


® r 

SduS 


PQ(u+l) 


PG(u-l) 


Er(k,n(k,u))  + £‘il  (k,n(k,u)) 


1 r2pou 


+ l du 


O JO 


± \ 

:,n(k,u)  + e11(k,n(k,u))  j ’ 


(32b) 


JUk.u)  = 2~(2pQku-k2) 


i. 


- We  present  plots  of  |Et  (p)/E  | vrs.  p/p  for  various  r (Graph  1). 

p-h 

Since  it  is  quite  difficult  to  obtain  a manageable  analytic  approximation  to  the 
integrals  in  Eq.  (28)  for  intermediate  momenta  (p  ~ PQ) » we  have  resorted 
where  necessary  to  numerical  integration.  From  the  numerical  integration  we 
have  shown  that  the  p4  dependence  for  small  p is  valid  up  to  ~.lp  . Increasing 
p further,  the  curves  rise  somewhat  rapidly  to  a relatively  broad  maximum  at 
p ~2pQ  (with  values  of  Ej  (p)/EQ  at  the  maximum  ranging  from  -.3  for  rs“3 


to  -.05  for  rg=2)  and  then  fall  rather  slowly,  droping  only  ~102  in  going 
from  p - 2po  to  p - 6pQ.  We  found  numerically  that  the  ^ variation  for 


p » pQ  becomes  discernable  for  p > - 20pQ. 


IV.  PLASMON  CONTRIBUTION 

* — * 

As  noted  earlier,  the  plasmon  contribution  to  Ej(p)  = Im  E (p)  arises  when 
the  pole  of  e 1 (k,E°(k)-E°(k-p)  +iri)  falls  within  the  domain  of  integration 
(Eq.  (27)). 

We  first  address  the  question  of  whether  a positron  of  momentum  p < pQ  is 

in  principle  able  to  excite  a (real)  plasmon.  Of  course  we  are  assuming  that 

2 

the  system  has  plasmons  of  energy  less  than  for  otherwise  excitation  would 

not  be  possible  on  energetic  grounds  alone  . From  the  geometric  construction 

of  Figure  (4a)  for  the  case  p < PQ,  it  is  apparent  that  Im  e(k,  E°(k)-E°(p-k)  +irj) 

* 

is  nonzero  throughout  the  region  of  integration  for  E^  (p).  Thus  the  positron 
necessarily  will  excite  only  electron-hole  pairs,  i.e.,  it  is  not  possible  for 
a positron  cf  momentum  p < pQ  to  excite  a plasmon,  regardless  of  how  small  the 
plasmon  energy  may  be. 


17 


We  turn  now  to  the  case  p > pQ.  From  Eqs.  (24),  (27)  and  the  comment 
after  (25),  we  see  that  the  poles  of  the  inverse  dielectric  function  are 
encountered  during  the  k-integration  for  E^*(p)  when  both 


and 


E°(p)  - E°(p-k)  = U>p(k) 


k < k 


(33a) 


(33b) 


i.e.,  when  it  is  dynamically  possible  for  a particle  to  excite  a plasmon. 

Let  us  for  the  time  being  neglect  the  small  dispersion  of  the  plasmon  mode 
and  set 


<i)  (k)  = a)  (o)  = a) 
p p o 


(34) 


whence  the  conditions  (33)  for  a pole  become: 


( P 2 ~ | P~k ) 2)  = CD 


2m 


k <k 


with  the  cutoff  momentum  in  the  absence  of  dispersion,  kc  given  by: 


(35a) 


(35b) 


k ° = (p2  + 2mo)  )Js  - p 

c ro  o *o 


(36) 


The  set  of  points  k satisfying  Eq.  (35a)  is  a sphere  centered  at  k ■ p 

and  of  radius  (p  - 2mo)o)  . When  (p2  - 2mo)0)^  < pQ,  points  k on  this  sphere 

are  such  that  k > kc°,  and  thus  do  not  satisfy  the  second  requirement  for  a 

pole,  Eq.  (35b).  | Alternatively  put,  as  is  evident  from  Figure  (4a),  these 

* 

points  lie  within  a region  of  the  domain  of  integration  for  Ej  (p)  in  which 


Im  e(k,  E°(p)-E°(p-k)  + in)  + 0 


Ej  (p)  in  this  case. 


There  then  is  no  plasmon  contribution  to 
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On  the  other  hand,  for  radius  (p2  - 2mu)oy  > pQ,  the  sphere  of  points 
satisfying  Eq.  -<35a)  does  at  least  pass  through  regions  of  the  domain  of  inte- 
gration in  which  Im  e(k,  E° (p)-E° (p-k)  +i6)  =*  0,  namely,  referring  to  Figure 
(4b),  the  crescent  shaped  region  bounded  by  both  the  outer  branch  of  the  Pascal 
limacon  and  the  sphere  of  radius  p,  and  the  oval  region  bounded  by  the  inner 
branch  of  the  limacon.  Only  in  the  oval  region,  however,  do  the  points  on  the 
sphere  also  satisfy  the  additional  condition  for  a pole,  Eq.  (35b) 

(Re  e(k,  E°(p)-E°(p-k)  +ir|)  0 in  the  crescent  region  ) and  it  is  therefore 

It 

from  only  this  region  that  we  get  a plasmon  contribution  for  Ej  (p). 

From  this  discussion,  we  conclude  that  the  value  of  positron  momentum  at 
the  threshold  for  excitation  of  a plasmon  is,  neglecting  dispersion,  given  by 

pth  = (2nWo  + Po^  * (37) 

2 

i.e.,  the  positron  has  at  threshold  an  energy  ^ such  that  upon  loosing  an 

P2 

energy  a)Q  to  the  plasmon,  it  ends  up  with  the  Fermi  energy  2^~ • 

Precisely  the  same  value  of  plasmon  excitation  threshold  momentum  is  found 
for  the  case  of  an  electron  probe  interacting  with  an  electron  gas  jagain 
neglecting  plasmon  dispersion  (Quinn,  1962 )j  . One  might  offhand  say  that 
the  value  p°^  is  reasonable  for  the  electron  case  since  the  external  electron 
cannot,  after  exciting  a plasmon,  end  up  with  an  energy  less  than  Eq,  the  corres- 
ponding states  already  being  occupied.  This  line  of  reasoning  would  seem  to  lead 
to  the  following  "paradox."  For  a positron  scattering  from  an  electron  gas,  one 
is  then  tempted  to  say  that  the  plasmon  excitation  threshold  should  occur  at  E°(p)  = 
u)Q,  with  the  scattered  positron,  being  distinguishable  from  the  electrons,  ending 
up  in  the  E°(p)  « 0 state. 

The  discrepancy  suggested  by  the  last  argument  lies  of  course  in  its  failing 
to  take  account  of  the  implications  of  momentum  conservation.  If  a positron  of 
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energy  E (p)  = ui  were  to  excite  a plasmon,  it  would  have  to  give  up  a momentum 


% o W 

(2110^)  . But  this  is  in  excess  of  the  plasmon  cutoff  momentum  k^  = (p*+2mu>o)^ 


- pQ,  therefore  implying  that  such  a positron  is  after  all  unable  to  excite  a 


plasmon. 


Note  in  fact  that  at  the  actual  threshold  pc^  (Eq.  (37))  the  plasmon  excited 


will  have  momentum  equal  to  the  cutoff  momentum  . As  positron  momentum  is 


increased  above  P^*  it  will  be  able  to  excite  plasmons  with  momenta  ranging 
from  a low  of  p - (p*-2mu)  )**  = k [forward  scattering]  to  a high  of  k ° |the 

p . 

positron  scattering  at  an  angle  cos  — from  its  original  direction  . These 

P i 


statements  follow  simply  from  Figures  (3c  & d) . Note  that  k^in  tends  to  zero 


as  p tends  to  infinity 
o 


For  p >p ^ we  thus  have  a plasmon  contribution  to  (p)  which  we  denote 


by  I,  (p),  and  which  is  given  by: 

V 


Et*  (p)  - - — 


p£ 


2tt 


/ 


d3k  Im 


v(k) 


oval 


e(k,  E°(p)-E°(p-k)  +in) 


th 


(38) 


In  the  oval  region  we  approximate  the  dielectric  function  by  its  high  frequency. 


low  wave  vactor  form  (i.e.,  again  neglecting  dispersion). 


e(k,u)  -1-22. 

u)2 


(39) 


Upon  using  Eq.  (39)  in  Eq.  (38),  introducing  spherical  coordinates  with  p taken 
along  the  polar  axis,  and  performing  the  trivial  azimutal  integration,  we  arrive 
at: 


px-pr 


*1*  (P)  = • tr Im  / dx  / 

rt  n 


dk 


k2(2px-k) 2 


(40) 


(-|  k(2px-k)  +1^/  -M2 


1 
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where  x ScosQ,  0 being  the  polar  angle. 

Doing  first  the  k-integrat'ion,  only  one  of  the  four  simple  poles  of  the 
integrand  falls  within  the  region  of  integration,  this  being  the  one  at 


px  - (p2x2  + 2mu)  - iri) 


(41) 


Applying  the  symbolic  identity  (for  real  x) 


STTTF  ■ p I + 


(42) 


it  is  trivially  seen  that  the  only  contribution  to  the  imaginary  part  of  the 
integral  arises  from  the  residue  at  this  pole.  Proceeding  in  a straight 
forward  manner,  we  eventually  find: 


VP)  = 


p2  + 2mcj  )**  - p 
o o ro 

p - (p  - 2mU)o)^' 


th 


(43) 

We  should  remark  that  precisely  the  same  result  is  obtained  for  the  plasmon 
contribution  to  the  imaginary  part  of  the  electron  self  energy  in  an  electron 
gas  (see  Quinn  (1962)).  Though  this  is  consistent  with  our  earlier  observation 
of  the  equality  of  the  positron  and  electron  plasmon  excitation  threshold  mo- 
menta, the  equality  of  the  two  self-energy  contributions  is  slightly  surprising. 

* 

We  note  the  following  about  the  behavior  of  ET  (p) . Recall  that  this 

pit  * 

expression  applies  only  when  p > p°  » p + k°  . Of  course  at  p = p°,  , IT  (p) 

cn  o c tn  I 

* o P* 

must  (and  does*)  vanish.  (p)  drops  abruptly  just  beyond  P,.^  with  a finite 

P£ 

derivative  ai.r.t.  p at  p » p°.+. 

tn 


* Note  that  the  numerator  in  the  argument  of  the  log  in  (43)  is  precisely  kc  . 
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•:;.’TT"EESSKS 


For  p » p , 
r ro  * 


(p)  ® - ~~  -2-  (r^-  "*■  00 ) 

V aoP  p0  po 


As  does  the  particle-hole  contribution  to  E^  (p) , the  plasmon  contribution 
then  decreases  with  increasing  positron  momentum. 


In  view  of 


Eqs.  (32)  and  (44),  and  using  E (p)  = £_  (p)  + E (p) , 

I I 1 o 


we  thus  arrive  at  the  sensible  conclusion  that  the  total  collision  rate  of  a 
positron  moving  in  an  electron  gas  decreases  with  increasing  momentum.  |i.e., 
that  for  increasingly  high  momentum,  the  particle  behaves  as  an  increasingly 
nondisturbing  probe  of  the  electron  gasj.  Moreover,  comparing  Eq.  (32)  and  Eq. 
(44)  we  see  that  the  plasmon  contribution  to  the  collision  rate  dominates  at 
high  momentum.  Thus, 


y(P)^2jEI*  (p)  | 2 j E*  (P)|- 

pi, 


-2.  £n  -2-  , -E- 

a0p  Po  po 


Let  us  remark  that  though  this  last  result  was  derived  for  a charged 

particle,  which  is  distinguishable  from  the  electrons  in  the  gas,  the  same 

result  applies  even  to  an  energetic,  "external"  electron,  since  at  very  high 

momentum  it  "loses  sight"  of  the  Pauli  restrictions.  These  restrictions 

would  only  come  into  play  if  the  electron  were  to  "try"  to  undergo  a highly 

unlikely,  strongly  inelastic  collision  which  would,  in  the  absence  of  the  Pauli 

restrictions,  have  brought  it  down  below  the  Fermi  sphere.  Indeed,  the  particle- 

hole  contributions  to  the  electron  and  positron  self-energies  should  thus  (and 

do)  approach  one  another  in  the  high  momentum  limit. 

* 

Before  presenting  numerical  results  for  E (p)  we  make  a worthwhile  modifi- 

V 

cation  of  the  formula  Eq.  (43).  Recall  that  we  have  until  now  neglected  plasmon 
dispersion.  This  was  done  for  reasons  of  calculational  expediency  and  was 

* 

It  can  be  shown  that  the  sign  of  the  particle  charge  does  not  affect  the  result, 
Eq.  (44). 
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justified  since  the  dispersion,  being  relatively  weak,  should  not  enter  into 
the  self  energy  in  a way  that  would  qualitatively  alter  any  conclusion  based 
on  neglecting  it.  However  there  is  one  quantitative  aspect  which  we  would  like 
to  be  (and  which  is  not)  contained  in  the  form  Eq.  (43)  for  E (p) , namely, 

V 

the  correct  value  of  particle  momentum  at  the  plasmon  excitation  threshold. 

We  have  seen  earlier  that  this  occurs  at  p = p°,  = p + k 0 if  the  dispersion 

is  neglected.  Now  allowing  for  dispersion,  at  the  plasmon  excitation  threshold, 

the  plasmon  produced  will  have  the  different  plasmon  cutoff  value,  kc,  for  its 

momentum.  However,  at  this  high  momentum  the  plasmon  energy  Wp(kc)  will  be 

somewhat  greater  than  01^(0)  = (usually  by  « 20%  for  most  metals).  Since  at 

* 

threshold  the  positron  still  ends  up  with  momentum  pQ  after  excitation  , it  must 
now  have  an  initial  momentum  greater  than  the  previous  value  P°jj. 

We  see  then  that  account  must  be  taken  of  dispersion  in  order  to  obtain  the 
correct  threshold  momentum.  On  the  other  hand,  for  particle  momentum  much 
higher  than  threshold,  the  plasmons  which  are  excited  (mostly  in  the  near  forward 
direction)  will  have  small  momenta,  and  energy  close  to 
E°(p)  - E°(p-k)  = t»)p(k)  is  satisfied  by  k + 0 as  p + since  tdp(k)  constant 


as 


k -+■  oj.  Thus  in  the  high  p regime,  the  result  Eq.  (43)  is  quantitatively 
correct  (within  our  overall  approximation). 

A device  for  shifting  the  threshold  momentum  to  the  correct  value  while 
preserving  the  correct  high  p behavior  of  the  formula  Eq.  (43)  is  to  everywhere 
replace  ti)Q  there  by  the  expression  Wp(k(p))  : 


E°(p)  - E°(p  - kc)  = wp(kc)  = (from  Eq.  (26))^  (k*  + kc  pQ)  f Pth  * PQ  + kc- 
th  th 
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* 


lT*  (p) 


u>  (k(p)) 


PS. 


2aoP 


Jin 


|(p2  + 2m  top(k(p))) 


p - Cp2  - 2m  u>  (k(p)))' 


k(p)  < k 


(46) 


Here  k(p)  (<k£)  satisfies 


E(p)  - E(p-k)  = a) p (k) 


(47) 


where  we  set  k | | p.  Note  that  the  important  k's  actually  are  parallel  to  p 


both  for  p near  p . and  for  p ■+•  °°. 
th 


I 

* 


It  is  seen  that  for  p at  the  proper  threshold  p , = k + p , E (p) 


'th  c ' I , 

* , p* 
vanishes,  as  desired  . Moreover,  since  k(p)  -*■  0 for  high  p,  we  obviously 


recover  in  this  range  the  earlier  correct  result  Eq.  (43).  The  expression 


Eq.  (46)  then  will  serve  as  our  approximation  to  E^  (p) . 


P£ 


In  practice  we  will  make  the  quadratic  approximation  for  ojp(k) 


3 Fo  ,.2 


OJ  (k)  — (jJ  + -Jr  — — k* 

p O 10  2 

v m a) 

O 


(48) 


Using  Eq.  (48)  in  Eq.  (47)  we  find  k(p)  which  is  then  substituted  back  into 


u)p(k)  to  give  C0p(p)  : 


,(k(p))  = ojo  + 8 |^P2-  (2mB+l)  2mo)o  - 2p  (p2- (2mB+l)  h | 


(49) 


with. 


B = 3 0 

B ' 10  2 


(50) 


m (i) 


In  this  approximation,  the  cutoff  momentum  kc  is  given  by: 


i*-(Fsy 


+ 4 “o  S - 8 


M 


(51) 


The  numerator  of  the  argument  of  the  logarithm  is  precisely  kc- 


' 
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Values  of  the  original  and  corrected  threshold  momentum  for  relevant  values 


of  r appear  in  Table  (1).  Differences  of  the  order  of  ~ "10Z  are  seen.  A 
s * 

comparison  for  r =■  2 of  E (p)  and  ET  (p)  fusing  Eq.  (49)  in  the  latter 
1 * 1 Pi  pi  1 J 

appears  in  Graph  (2).  It  is  apparent  that  both  results  are  even  quantitatively 

very  similar,  save  for  the  discrepancy  in  thresholds.  We  point  out  that  deriva- 

i 

* 

tive  oj.r.t.  p of  ET  (p)  is  again  nonzero  at  p * p , +. 

l . • tn 

pi  i 

Finally,  we  present  plots  of  |Ej  (p) | /E  vrs  p/pQ  for  relevant  r 

pH 

Graph  (3)  • The  following  qualitative  features  emerge:  an  abrupt  rise  in 
* 

|Ej  (p)|/EQ  just  beyond  p = pt>  , a fairly  broad  maximum  of  magnitude  ~ 1/10  at 
Pi 

p/pQ  - 2-4,  followed  by  a rather  gradual  decrease  to  zero. 

it  it  it 

The  total  value  for  E (p)  is  given  by  the  sum  ET  (p)  + Ey  (p) 

1 p-h  pf. 

(Eqs.  (29)  and  (46)).  Before  proceeding  to  a discussion  of  the  character  of 
* 

Ej(p)  it  is  necessary  that  we  verify  the  notion  implicit  until  this  point 

that  the  probe  interacting  with  the  electron  gas  forms  a reasonably  well  defined 

quasiparticle  state.  More  precisely  put,  we  must  verify  that  the  energy  width 

* 

of  our  quasi-eigenstates,  namely  E^(p)  is  small  compared  to  the  excitation 

energy  E°(p).  Plots  of  E ^p/p^) /E°  (p/p^)  vrs  p/p^  are  presented  in  Graph  (4) 

for  representative  r . It  is  seen  that  only  rarely  does  this  ratio  approach 

s 

even  a tolerable  - 20Z.  For  the  most  part,  the  fractional  uncertainty  is 
considerably  smaller.  We  can  be  assured  then  that  our  results  are  reasonably 
sensible. 

We  are  now  able  to  present  the  following  picture  for  the  RPA  total  collision 
rate  of  a positron  moving  through  an  electron  gas.  |see  Eqs.  (29)  and  (46) j . 
Plots  of  |Ej(p)|/EQ  * y/Eq  vrs.  p/pQ  appear  in  Graph  (5)  for  relevant  rg. 

For  p < pQ,  only  particle-hole  pairs  are  excited.  |Ej(p)]  starts  off 
with  a - Pl*r£.2  dependence  and  increases  smoothly  through  the  region  p = pQ. 


At  a momentum  p = p . - 2p  (for  r of  interest)  the  particle  begins  to  excite 
tn  os 

plasmons  |in  addition  to  p-h  pairs J , this  being  manifested  in  a definite  jump 

in  |E^(p)|  there.  There  is  then  a rather  broad  maximum,  of  magnitude  ranging 

from  - .2  E for  r = 1 to  ■ .7  E for  r =3,  reached  at  p ' 2-4  p . This 
os  os’  r *o 

is  followed  by  a relatively  slow  drop  off  with  an  asymptotic  ~ — £n 

P ^o 

dependence,  originating  from  the  dominating  plasmon  excitation  channel. 

For  reference  and  convenience  we  present  plots  of  E^  (E°)  vrs  E°  in  Graph  (6). 

★ 

In  Graph  (7)  we  present  a comparison  between  our  results  for  E^(p)  for  a 
positron  and  Quinn's  results  for  the  same  for  an  electron.  The  larger  magnitude 
of  the  imaginary  part  of  the  positron  self  energy  is  readily  understood: 

Since  the  positron  does  not  see  the  Pauli  restrictions  operative  for  the  electron, 
it  has  more  available  phase  space  to  scatter  into, resulting  in  its  having  a higher 
collision  rate. 

The  mean  free  path  A(p)  of  a positron  quasi-particle  with  excitation  energy 


E(p)  may  be  taken  as 


A(P)  -iS&l 


Now  E(p)  is  given  approximately  by  Eq.  (11a).  We  will  here  ingnore  the 
k-dependence  of  Re  E (p,E(p))  and  take  simply: 

‘ > ' TTFT  1 

A plot  of  A(E°)  vrs  E°  appears  in  Graph  (8).  The  correspondence  between  the 


structure  in  Ej  described  above  and  that  in  A is  apparent. 


For  E (p)  « E . 
r o 


A (E°) 


/2  15 


2 3/2l  3/2 

e2ao7T  » / Eo 


I 


1 


(55) 


and  for  E (p)  » E 


^E°>  “ 4ao  J- 


o 2,n(E°/E  ) in  E° 
o 


A minimum  of  the  order  of  a few  Angstrom's  occurs  at  E *•  2Eq. 

The  author  would  like  to  thank  Professor  Walter  Kohn  for  numerous  valuable 
discussions. 
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TABLE  1 


r 

s 

Po 

E 

o 

k o 
c 

k 

c 

vv 

o 

pth 

pth 

(a.u. ) 

(eV) 

po 

po 

(eV) 

(eV) 

m 

Po 

2 

.960 

.12.5 

.528 

.640 

16.7 

21.4 

1.53 

1.64 

3 

.640 

5.57 

.621 

.739 

9.07 

11.3 

1.62 

1.74 

4 

.480 

3.13 

.698 

.820 

5.89 

7.23 

1.70 

1.82 

5 

.384 

2.00 

.764 

.882 

4.22 

5.11 

1.76 

1.88 

a = Bohr  radius 
o 


t»2 


me 


2 


v 


Pq  = Fermi  momentum  =» 
Eq  = Fermi  energy  = 


(atomic  units) 


k ° = plasmon  cutoff  momentum,  neglecting  dispersion  = (p1*  + 2irto  )**  - p 


k£  = plasmon  cutoff  momentum,  including  quadratic  dispersion. 


fe  - (feV  - *.(&  -»)H 


3 Po 


10  2 
m a) 


U)p(o)  * plasmon  energy  at  cutoff,  neglecting  dispersion 


Wp(kc)  * plasmon  energy  at  cutoff,  including  quadratic  dispersion. 

” positron  momentum  at  plasmon  excitation  threshold,  neglecting  dispersion: 
p_.  = k + p . 

*th  c *o 

Ptjj  = positron  momentum  at  plasmon  excitation  threshold,  including  quadratic 
dispersion:  pt^  * kc  + p^. 
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FIGURE  1:  Contours  for  Self  Energy  Evaluation 


a)  Original  contour  C for  the 
frequency  integration, 
x denotes  poles  of  dielectric 
function. 

+ denotes  a pole  of  the  free 
positron  Green's  function. 

+ occurs  in  the  1st  quadrant 
when  a)  > E°(p-k),  e.g. , + 
and  it  occurs  in  the  2nd 
quadrant  when  E°(p-k)  > u>» 
e.g.,  +2- 

I 

j 

b)  Close  contour  in,  e.g., 
upper  half  plane. 


c)  Rotate  contour  in  b 
by  + 90°.  A residue 
contribution  arises  when 
the  pole  + lies  in  the 
1st  quadrant.  The  contour 
C'  may  be  now  deformed  into 
the  imaginary  axis. 


+2 


FIGURE  3:  Cross  Section  of  Region  of  k-Space  Contributing  to2->T  (p) 

ph 

Figure  a,  p < pq;  Figure  b,  p > pQ.  0 is  the  origin  of  k-space.  Shaded 

regions  correspond  to  points  k lying  both  within  the  domain  of  integration 

in  Eq.  (21)  (within  the  sphere  of  radius  p centered  at  k = p)  and  within 

1 2 

the  region  of  k-space  in  which  Im  e(k,-^(.2  pk  cos  9 - k ) 
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Plasmon-Contribution  to 


Im^Z  (p) 


dotted  spheres  in  both  cases.  Only  when  its  radius  is  greater  than  pQ  do  the  points 
on  this  sphere  lie  in  regions  where  Im  e(k,  (2  kp  cos  0 - k ))  =0  (b),  and  then, 
only  points  such  as  kj  in  the  small  oval  give  a zero  of  £(k,  (2  kp  cos  9 - k2)). 

(Re  £(k,  (2  kp  cos  9 - kz))  ^ 0 for  points  such  as  k ) . Obviously,  threshold 

corresponds  to  the  p for  which  the  dotted  sphere  coincides  with  the  sphere  centered 
at  k = p and  of  radius  p . 
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